Number-Theoretic Nature of Communication in Quantum Spin Systems 
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Abstract 

The last decade has witnessed substantial interest in protocols for transferring information on networks of quantum mechanical 
objects. A variety of control methods and network topologies have been proposed, on the basis that transfer with perfect fidelity 
— i.e. deterministic and without information loss — is impossible through unmodulated spin chains with more than a few 
particles. Solving the original problem formulated by Bose [Phys. Rev. Lett. 91, 207901 (2003)], we determine the exact 
number of qubits in unmodulated chains (with XY Hamiltonian) that permit the transfer with fidelity arbitrarily close to 1, a 
phenomenon called pretty good state transfer. We prove that this happens if and only if the number of nodes is n = p— 1, 2p — 1, 
where p is a prime, or n = 2 m — 1. The result highlights the potential of quantum spin system dynamics for reinterpreting 
questions about the arithmetic structure of integers, and, in this case, primality. 
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Introduction. - Since the pioneering work by Bose 
Q , quantum channels implemented by spin systems have 
been interpreted as wires for the transmission of possi- 
bly unknown qubit states. Motivated by the perspective 
of designing quantum buses and nanodevices, a vigorous 
effort has delineated the field of quantum spin systems 
engineering In this context, one of the desirable tasks 
is to transfer the state of a particle into another one with 
maximum fidelity; when the fidelity is 1, we have perfect 
state transfer (PST), a notion originally introduced in 
0- Given an inherent difficulty in the manipulation of 
coherent quantum mechanical objects, the most appeal- 
ing set up for state transfer employs a time-independent 
Hamiltonian and no interaction with the system except 
at initialisation and read-out. If we do not use ad hoc 



coupling schemes and control protocols (see [14j ; also the 
recent [23| and the references therein), a spin chain (ID 
magnet), with a Heiscnbcrg XY Hamiltonian, exhibits 
PST only between 2 and 3 qubits 0- This is a nega- 
tive result because such a chain is arguably the quantum 
wire obtainable with the smallest amount of physical and 
technological resources. 

On the basis of sufficient conditions for PST, this fact 
fueled a mathematically challenging classification pro- 
gramme aimed to identify PST in general network topolo- 
gies [ll|. The £;-dimensional hypercube (with n = 2 k 
qubits) has the best known performance in terms of the 
distance travelled by a single excitation: PST occurs be- 
tween antipodal nodes at (network) distance k. However, 
for practical purposes, it is natural to study whether a 
minimalist structure like the chain (with n qubits) can 
still be usefully employed. In particular, even if we al- 
ready know that there is no PST for n > 4 we may 
still ask the following question: given n and an e > 0, 
does there exist t such that the fidelity at time t between 
qubits 1 and n is larger than 1 — el When the answer 



is "yes" , we say that there is pretty good state transfer 
(PGST). While Bose 0] verified that the fidelity could 
be remarkably high even for rather long chains, the no- 
tion of PGST was formally isolated in [10( as a relaxation 
of PST. We give, in the present Letter, a complete char- 
acterization of the parameters for which there is PGST. 
Our findings can be combined into the following simple- 
sounding statement: 

Theorem. - A uniformly coupled chain of n particles with 
XY Hamiltonian has PGST if and only if n = p — 1 or 
2p — 1, where p is a prime, or if n — 2 m — 1. 

The significance of the result is twofold. From the 
physical point of view, it is valuable that we rigorously 
describe a phenomenon with applications to the study of 
quantum nanodevices Q. In fact, the related observa- 
tions obtained up to now are either numerical or frag- 
mented. The message conveyed by the notion of PGST 
is that an unknown qubit state can be transferred with 
arbitrarily large fidelity between the end nodes of "long" 
chains, whenever the waiting time is not an issue. From 
the mathematical point of view, we highlight properties 
of quantum interference as a consequence of number- 
theoretic constraints. The appearance of prime numbers 
in the theorem indicates a connection between the occur- 
rence of PGST in chains of a given length and primality 
testing. This suggests a potential application of quantum 
dynamics on graphs for reinterpreting number-theoretic 
problems, when the problem description is encoded in 
the parameters of the system. 

It is well established that the dynamics of a single exci- 
tation in networks of spins with unmodulated couplings 
is a continuous-time quantum walk on the unweighted 
graph modeling the network. Hence, our result can also 
be reinterpreted in the language of continuous-time quan- 
tum walks: it describes the maximum hitting probability 
between antipodal vertices induced by a one-dimensional 
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quantum walk. Walks of this type have been studied 
in great detail [161 ]. This analogy prompts us to disre- 
gard the original spin system setup and to work with 
a single n-level system. The realisation and simulation 
of this quantum device is the centre of several discus- 



sions [21] . It is remarkable that at a wider level, state 
transfer can be seen as the simplest model in a family of 
processes for quantum transport: for example, the pop- 
ulation transfer of n- level systems [22|], after including 
the role of the environment, the transport of energy in 
organic molecules [2(j, and a closely related mathemati- 
cal setting describes state transfer in chains of harmonic 
oscillators coupled with beam-splitter interactions flij . 
Finally, the experimental implementation of the Heisen- 
berg XY chain has been proposed in schemes including 
cold-atom optical lattices and superconducting circuits 

1- 

PGST. - The Hamiltonian governing the evolution of 
the system acts on the Hilbert space H = Cf ® • • • ® 
C 2 . If we do not include external static potentials, the 
Hamiltonian is 



_^ 71 — 1 

u=l 

where er™ (w G {x, y, z}) is a Pauli matrix on C 2 and 
J u is the coupling strength between the particles u and 
U + 1. By virtue of the Jordan- Wigner transform 17j . 
the free evolution for a time f of a single excitation 
originally located at site \u) € {|1), |n)} is given by 
e lHt \u) = U(t)\u), where H UyV = J u if v = u + 1 or 
u = v + 1, and H u v — 0, otherwise. The n x n real 
symmetric matrix H is the Hamiltonian restricted to the 
single excitation sector. PST occurs between 1 and n 
if there is a t e K+ such that \(n\U{t)\l)\ = 1, i.e., 
the channel has maximum fidelity, PGST occurs between 
1 and n if for every e > there is t € IR + such that 
|(n|J/(t)|l)| > 1 — e. Analytical solutions for coupling de- 
sign able to achieve PST for any n have been presented 
in previous works (see [3]). Abstractly, a chain of length 
n is modeled by a network called an n-path and denoted 
by P n . The links representing the particle-particle inter- 
actions are {1, 2}, {2, 3}, {n — 1, n}. When the chosen 
couplings are uniform (w.l.o.g., J u — 1), the Hamiltonian 
H is the adjacency matrix of the n-path. (Recall that the 
adjacency matrix of a graph has ij-th entry 1 if there is a 
link between the nodes i and j; 0, otherwise.) We know 
from 10] that there is PGST in P 4 and P 5 . 

We shall first prove the theorem. Next, we present 
some details about the cases when there is no PGST. In 
particular, we will give an explicit upper bound on the 
fidelity in a special case. In Appendix 3, we discuss PGST 
between internal nodes, by considering a link between 
state transfer and control theory on networks (see 0]). 

Proof of the theorem. — The proof of the theorem is 
based on a direct linear-algebraic analysis of the eigensys- 



tem of U (t) and on the application of standard number- 
theoretic tools, especially Kronecker's theorem on Dio- 
phantine approximation. We begin by considering a gen- 
eral property of PGST and basic facts about bipartite 
graphs. We then use relations on the eigenvalues with a 
detailed case- by-case treatment. 

For any two vertices u and v of a graph, 

U{t)\u) ~ 7 |«> = 7 C/(<)(7" 1 |w) - U{-t)\v)); 

here 7~ 1 /— U{— t) is the Hermitian adjoint of r )^ l I—U{t). 
Since U(t) is unitary and ||7|| = 1, we see that 

\\U(t)\u)- 7 \v)\\ = \\U(t)\v)- 7 \u)\\. 

So, if we have PGST from u to v, we also have it from v 
to u. 

A graph X is bipartite if there is a bipartition of the 
set of nodes such that the links connect only nodes in 
different parts. Suppose X is bipartite and let D be a 
diagonal matrix such that D u u is 1 or — 1, accordingy 
as u is in one or the other part of the bipartition. Then 
DAD = - A, and if U(t)\u) w 7 |w), 

jD\v) « DU{t)DD\u) = U{-t)D\u). 

But \u) and \v) are eigenvectors for D with eigenvalues 
1 or —1; the eigenvalues are equal if and only if u and v 
are in the same part. So there is a sign factor a u>v and 
l\v) « o- u . v U(-t)\u). Accordingly, U(t)\v) « 7~ 1 er„ il ,|u). 

By the above, U(t)\v) ~ 7|i>), and we conclude that 
7 sa 7~ 1 er u .„. Hence, 7 ~ ±1 if u and v are in the same 
part, and 7 w ±i if they are not. (For PST this obser- 
vation is due to Kay Q.) Let F denote the permutation 
matrix of order n x n such that F\r) = \n + 1 — r) for 
all r. Let E\, . . . , E n be the idempotents in the spectral 
decomposition of the path P n [i.e., of its adjacency ma- 
trix). We can then write F — Y^-=i{~^Y~ 1 Er- If we 
have PGST at time t then U(t) « 7F and, therefore, 

1 = detE/"(t) «7 n det(F) = 7 ™(-l)L"/ 2 J . 

This yields three cases: (1) if n = 1 mod 4 then 
(_1)L«/2J = 1 and 7 w 1; (2) if n = 3 mod 4 then 
(_1)L«/2J = _i anc j 7 ~ _ 1; (fy if n is even t hen 

i n = (— l)"/ 2 and 7 ~ ±i. It is well known that the 
eigenvalues 8 r of P n are given by 6 r = 2cos[7rr/(n + 1)]. 

We start with the positive results. If PGST occurs then 
U(t) gets arbitrarily close to 7F. This means that e l0rt ~ 
( — l) r_1 7 for r = 1, ...,n. Set m = [n/2\. Assume 
7 = ±1 if n is odd and ±i if n is even. 

First, we prove that for the path P n , if e 16 ^ ~ 
(-l) r - 1 7, for r = l,...,m, then e iBrt « (-l) r ~ 1 7, for all 
r = 1, and, hence, C/(t) w 7F. To see this, assume 
e i9 T t py ^_j^r-i^_ Since 71-paths are bipartite, 0„ + i_ r = 

For 



6* r , and, therefore, e* 1 



(-irv 



PGST, we need (-l) n - r 7 = (-l) 7 -^" 1 , or, equiva- 
lently, 7 2 = (— 1)" _1 . As this holds for our choice of 
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7, we are done. Since 9 r = 2cos[7rr/(n + 1)] we have that 
9 r is a polynomial with rational coefficients of degree r 
in 6\. We also set 9q — 2. It follows that the numbers 
9q, . . . ,9<i are linearly independent over Q if and only if 
the degree of the algebraic integer 9\ is greater than d. 

This provides us with the necessary tools to prove the 
first part of the theorem: if n = p — 1 or 2p — 1, where 
p is prime, or if n = 2 m — 1 then we have PGST on P n . 
If a E C we use use Q(a) to denote the field obtained 
by adjoining a to Q. If F is a subfield of Q(a), then 
Q(a) is a vector space over F; its dimension is the index 
of F in Q(a). (For details see, e.g., [13], Chap. 5.) Let 
9 = 2cos[7r/(n + l)] and set C = e i7r /0 +1 ). Then 9 G Q(C) 
and £ is a root of the quadratic x 2 — x9 + 1. So the index 
of Q(0) in Q(C) is at most two. If n > 3 though, 9 is 
real and C is not. Thus, the index is exactly two. The 
degree of £ is 0(2n + 2), where 4> is Euler's function, and 
therefore, the degree of 9 is <f)(2n + 2)/2. If n = p — 1, 
2p—l, or 2 m — 1, where p is a prime, we have, respectively 

0(2n + 2) = 0(2p)=0(p)=p-l, 
0(2n + 2)/2 = 0(4p)/2 - = p - 1, 
<?(>[2(n + l)]/2 = 2 m - 1 . 

Finally, in each of these three cases the positive eigen- 
values of P n are linearly independent over Q. Next, let 
8 equal 0, |, j, or |, accordingly 7 equals 1, — 1, i, —i, re- 
spectively. For each r = 1, . . . , m, let a r = \ if r is 
even and a r = if r is odd. By Kronecker's theorem 
(see [HI), for any e,T > there is a t c > T, and inte- 
gers p r , r = 1, . . . , m, such that |i e 4f — p r — a r — 5| < e. 
It follows that as e -> 0, e^*'^ -> (-l) 1 *"^, with 
r = 1, . . . , n, so that we have PGST. Let us now focus 
on the negative cases: if n = mp — 1, where p is odd 
and m > 3, then PGST does not occur on P n . Suppose 
n + 1 = mp, where p is odd. After some algebra (see 
Appendix 1), we obtain the following equation: 

2 

{9 1 — 9 2 ) + 1)' (6> mr +i — #rw+2) 

r=l 

p-i 
2 

+ l) r (^mr-l — = 0. 

Denote the three terms on the left by Z), E, and F, re- 
spectively. If we have PGST, then there is a sequence 
of times (tk)k>o such that e l6rtk —> (— l) r_1 7, and so 
e i{8 s -e s+1 )t k _^ _ L T h ere f orei e i£>t* ^ _ 1) w hile e 4£; * fc 

and e lBtfc both tend to 1 or to -1. Thus, e^+s+^K- 
— 1, which is impossible, since + E 1 + F = 0. It is not 
difficult to verify that the cases considered above include 
all integers. This ends the proof of the theorem. A corol- 
lary is that if 2 < n < 10, we have always PGST, except 
in the case n = 8. 

Fig. Q] gives the smallest times needed to achieve a 
relatively large fidelity (> 0.99 and, thus, e = 0.01) for 
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FIG. 1: Logarithm of the earliest points in which the fidelity 
is strictly greater than 0.99 for 2 < n < 7. The numbers have 
been obtained by plotting |(7(t)i. n | in — 2, 7) and then by 
analysing sections of the curves. Clearly 0.99 is an arbitrary 
choice. Notice the jumps between the pairs (2, 3), (4, 5), (6, 7). 
Because PGST depends on the positive eigenvalues of P n , 
this phenomenon may be explained by the fact that \n/2\ = 
L(n + 1)/2J, for n even. 



chains of length 2 < n < 7. The parameters to be consid- 
ered for numerics are n and e. Indeed, the waiting time 
depends also on the tolerance e being close to one. The 
log plot suggests that for a fixed e there is a behaviour 
that is linear in n. 

It is a corollary of the theorem that there is PGST on 
P n if and only if its positive eigenvalues are linearly inde- 
pendent over the rationals. A proof of some cases when 
there is PGST can be constructed with the use of facts 
about linear independence of roots of unity discussed by 
Conway and Jones (in particular Theorems 1 and 7) and 
by Watkins and Zeitlin 0]. The crucial observation for 
PGST is indeed the linear independence (over the ratio- 
nals) of the numbers cos[7rj'/n + 1] for certain choices of 
3- 

Bounding fidelity. - We have seen that there are 
many cases in which there is no PGST for P n . We out- 
line now a general technique for proving upper bounds on 
the fidelity whenever this happens. Algebraic graph the- 
ory is again the natural toolbox to employ. The spectral 
decomposition of the adjacency matrix of a graph X is 
A = "^2 r 9 r E r . Two vertices u and v of X are cospectral 
if, for each r, the projections E r \u) and E r \v) have the 
same length. We say they are strongly cospectral if, for 
each r, we have E r \u) = ±E r \v). In [lfj, it is shown that 
if we have PGST (or PST) from u to v then u and v are 
strongly cospectral. If the eigenvalues of A are simple, 
two vertices are strongly cospectral if and only if they are 
cospectral. Assume U{t) u . v — ^ r {E r ) u . v e l9rt . Define e r 
by the requirement that (E r ) UjV — e r {E r ) u>u . (For paths, 
e r = (— l) r_1 .) Then U(t) UtV is a convex combination of 
the norm one complex numbers e r e lt ® T . For PGST to oc- 
cur, these numbers must all be approximately equal. We 
can see this by applying the triangle inequality, 

\U(t) 

u.v u.u\ = L 

r r 

In particular, PGST cannot happen if there is some set 
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S of eigenvalue indices such that 



u.u 



res 



A9 r t 



is bounded away from zero (for all t). With a direct anal- 
ysis of this expression, we can rule out PGST explicitly 
when n = 3k + 2 and k is even. (The details arc in Ap- 
pendix 2.) We leave open the challenge of finding explicit 
bounds on the fidelity in the remaining cases. The proof 
technique outlined here can be potentially extended to 
other network topologies. 

Conclusions. - By solving an open problem about 
quantum transport [3[, we have highlighted number- 
theoretic properties of quantum communication in spin 
chains. We have studied general properties of PGST. We 
have given necessary and sufficient mathematical condi- 
tions for PGST to occur on XY spin chains with uniform 
couplings. 

The physical intuition paralleling the mathematical re- 
sult suggests that the spin wave can reach the end of the 
chain with an arbitrary high peak only when the number 
of particles does not permit significant constructive in- 
terference. In this case, the trajectory of the amplitude 
(in a fixed time interval) suggests an intriguing analogy 
with chaotic dynamics that remains to be explored. 

Deciding whether there is PGST is computationally 
equivalent to primality testing, a task that is performed 
efficiently with the AKS test |l[ . Experimental detection 
of PGST would correspond to a natural algorithm for pri- 
mality. Its complexity would be determined via bounds 
on the time required by physical evolution (or its simu- 
lation) and on resources for tomography and sequential 
measurements. 

We have outlined a technique for quantifying maxi- 
mum fidelity when there is no PGST. Some cases remain 
with no complete answer. We have shown how notions 
of network control theory can be applied to study com- 
munication in spin systems. Exploring PGST in general 
networks beyond the n-path requires similar methods, 
but it is a challenging task. It is valuable to observe that 
we have considered a system without spatial disorder; its 
behaviour does not exhibit effects due to Anderson local- 
ization. Numerics in [j| indicated that the fidelity of this 
system tends to be robust when a relatively small amount 
of disorder is introduced in the couplings. On the other 
side, it was shown in [l5| that the speed of propagation of 
coherent walks is suppressed exponentially in the amount 
of imperfection. 

We have left open the development of a comprehensive 
theory of PGST. Such a theory is important to obtain a 
fuller understanding of transport in networks of quantum 
mechanical particles, either engineered or found in na- 
ture. Experimental tests based on photonic waveguides 
are currently investigated. (See [18j for background on 
such schemes.) 
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Appendix 1 - In the proof of the theorem, suppose 
n + 1 = mp, where p is odd. Then 

(p-i)/2 

1 + 2 £ (-l) r cos(7rr/p) = 0. 

r=l 

If we multiply this by cos (7r/ (n + 1)), we get 



cos 



j +B-D 



■nimr + 1)\ / ir(mr — 1) 

cos I — — ^ I + COS ' 



n+1 / \ n + 1 

which yields the following relations on eigenvalues: 



= 0, 



r=l 



r=l 



with a; = 1,2. The equation with a; = 2 is obtained if we 
multiply by cos (2ir/(n + 1)) the equation with x = 1. If 
we subtract these equations, we have 



(#1 — #2) + y^(-l) r (flmr+l - 

r=l 

p-i 
2 

+ ^(-l) r (fmr-l - C-2) = 0. 
r=l 



Appendix 2 - We have outlined a general technique for 
bounding fidelity. We know that there is no PGST when 
n = 3fc + 2 and fc is even. This is the case considered 
here. There is no PGST if there exists some set S of 
eigenvalue indices such that the expression 



£(£r) 



is bounded away from zero (for all t). In this case, Q\ — 
Ok + 8k+2- Define 



i8 r t 



Here, 



h(t) = X^Lne' 



(^r)l,n — (-l) r (^r)l,li 



where (-E r )i,i > and J2 r (E r )i,i = 1- So, we can rewrite 
in the form 



fe(t) = ^ flr (-irv 



Since ~^2 r a r = 1, if ~ 1 then the summands in this 
expression must be approximately equal. Consider the 
sum 

ai e ieit + (-l) fe -Ve i<M 
+ (-l) fe+1 a fe+2 e^+ 2t + (-l)" a „ e ^ 4 . 

For PGST to occur, its absolute value must be close to 
«i + flfe + «fc+2 + a„. We can simplify a little by working 
with 

+ (-l) fe+1 a fc+2 e j(efe + 2 - ei) * + (-l)"a„e 4(e "- ei) *; 

which has the same absolute value. We aim to show that 
the real part of this sum is bounded away from a\ + a k + 
+ a n . We note that 6 n — — #1, a n — a\, and each of 
a\, a k , afe+2, a„ is positive. 

When fc is even, we can write the real part of this sum 

as 

ai - a fe cos[(6» fe - 6>i)i] 

- a k+2 cos[(9 k+2 - Bi)t] - ai cos(2#ii). 

Since 9\ = 9 k + 9 k+2 , this is equal to 
ai - a k cos(9 k+2 t) 

- a k+2 cos(9 k t) - a x cos[2(# fc + 9 k+2 )t}. 

When 

cos(8k+2t),cos(0kt) < -v / 3/2, 

we have 

-1/2 < sin(0 fe+2 t),sin(6» fe t) < 1/2. 

Whence 

cos[(0 fe + 6 k+2 )t] > (3/4 - 1/4) = 1/2 

and so 

cos[2(6» fe + 9 k+2 )t] > -1/2. 

Consequently, 

ai - a k cos(8 k+2 t) - a k+2 cos(9 k t) 

- 01 cos[2(0 fe + 6 k+2 )t] 

< 3ai/2 + (a k + a fe+2 )v / 3/2 

= 2oi + a fe + a fe+2 - (ai + (2 - \/3)ak 

+ (2 - V3)a fe+2 )/2. 

On the other hand, if cos(9 k+2 t) > V3/2 then 

ai - a k cos(9 k+2 t) - a k+2 cos(9 k t) 
-ai cos[2(9 k + 9 k+2 )t] 
< 2a\ + V3a,k/2 + a k+2 + a\ 
= 2a\ + a k + a k+2 
- (2 - v / 3)a fe /2. 
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It follows that if cos(9kt) > V3/2 then we have an explicit 
upper bound: 

2ai +a k + a k+2 - (2 - Vs)a k+2 /2. 

This rules out PGST when n — 3k + 2 and k is even. 

Appendix 3 - Through this work we have studied 
PGST between the extremities of a chain. For practi- 
cal purposes it may be useful to have PGST between 
particles corresponding to internal nodes. We discuss an 
argument for showing that if there is PGST between in- 
ternal nodes then there is PGST between the extremities. 
In doing so, we appeal to algebraic techniques from quan- 
tum control theory of spin systems (Ref. [6]). Given a 
graph X with set of vertices V and adjacency matrix A, 
let z be the characteristic vector of some set S QV. We 
define and denote by 

W z = [zlAz]...^-^} 

an n x n matrix with entries in Z-°. The matrix W^is 
called the walk matrix of X with respect to S. The pair 
(X, z) is said to be controllable if the matrix W z is invert- 
ible (i.e., det(W z ) ^ 0). The set-up is a graph-theoretic 
analogue of the famous Kalman rank condition in control 
theory: the matrices e lAs and e lzz * (s,( e M + ) from a 
controllable pair generate a dense subgroup of the uni- 
tary group U (n) (n> 2). 

A closed walk is a sequence of vertices, with consecu- 
tive vertices adjacent, that starts and ends at the same 
vertex. If we have PGST on P n from k to £ then the 
vertices k and I are cospectral, so then k + £ = n + 1 and 
the generating functions for closed walks at k and at £ 



are equal. Consequently, the distance from k to an end 
vertex equals the distance from £ to the antipodal end 
vertex. 

By the A-module generated by a vector, we mean the 
smallest A-invariant subspace that contains the vector. 
If \£) lies in the A-module generated by and we have 
PGST from k to n + 1 - k, then we have PGST from £ 
ton+l-£. 

Suppose 

U(t)\k) RJ 7 |n+l - k). 

Because \£) lies in the A-module generated by there 
is a polynomial / such that 

\£) = f(A)\k). 

Then 

U{t)\£) = U{t)f{A)\k) - I(A)U(t)\k) « 7 /(A)e„ +1 _ fe . 

If $ is the 'flip' automorphism on the path, then $ 
commutes with A and so 

f(A)\n+l-k)=f(A)$\k) = 

$f(A)\k) =$\£) = \n + l-£). 

If n + 1 is a prime, then all vertices i give a controllable 
pair (P n ,i), and so PGST between any pair of vertices 
implies PGST between end vertices. If n + 1 = 2p where 
p is prime and k ^ 2,p, then PGST from k to n + 1 — k 
implies PGST between the end vertices. If k = p then k 
is the central vertex and PGST cannot occur at k. We 
leave open the case n = 2 m + 1. 
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